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Abstract. We compare two formalism that describe minijet production in 
pA and AA collisions: pQCD supplemented by Glauber-Gribov multiple semi- 
hard parton scatterings (pQCD+Glauber), and the Colour Glass Condensate 
(CGC). We argue that in a suitable limit they are equivalent to each other, the 
PQCD-|-Glauber model being more accurate from a numerical point of view. 
Finally, we analyze RHIC data on Au-Au integrated charged multiplicities in 
the pQCD-l-Glauber framework, and conclude that at least at central rapidity 
there is no sign of giuon saturation. 

Keywords: pQCD, multiparton scatterings, gluon saturation, CGC 
PACS: ll.80.La, 24.85. -fp, 25.75.-q 



1. pQCD, Glauber-Gribov parton rescatterings 
and the colour dipole 

The model of Ref. [1], hereafter labeled "pQCD+Glauber" , assumes minijet produc- 
tion in pA collisions to be dominated by semihard parton multiple scatterings. By 
"semihard scattering" we mean a process with a minimum transverse momentum 
exchange p^ ~ 1 ~ 2 GeV, described by leading order (LO) pQCD parton-parton 
cross section. The model assumes the S-matrix for a collision of one parton on n 
partons from the target to be factorizable in terms of S-matrices for parton-parton 
elastic-scattering, and assumes generalized pQCD factorization. Considering only 
gluons for simplicity, the minijet transverse spectrum is then given by 

J mj J A ^ P 

j^ = G{.yr)^ + AG{x'y^)^ , (1) 
aija pt a pt d'^pT 



where G{x, Q^) is the proton distribution function of a gluon with fractional momen- 
tum X, and x{x') « (j)t I \fs) exp(±77). In the second term of Eq. (1), the A-nucleus 
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partons are assumed to undergo a single scattering on the proton with cross section 

(2) 



The LO pQCD cross section for gluon-gluon scattering is -^p^ « i*^s (^^rp^^ ^^'^ 
the Hmits of integration on x[ are given by parton-parton kinematic constraints. 

The first term of Eq. (1) accounts for multiple semihard scatterings of the proton 
partons on the nucleus. The proton is considered pointlike at an impact parameter 
b. Nuclear effects arc assumed to be due only to multiple semihard scatterings, and 
are included in da'^^-'^'^/cPpT, the transverse momentum distribution of a proton 
parton who suffered at least one semihard scattering. This is written as [1-3]: 



n! J d-'ki d-^kn 

x5(2)(^k^_p^) , (3) 



n=l 



where (7p'"''^(po) = / <Pk Ji^'^'' is the integrated gluon-nucleon cross section, which 
depends explicitly on the infrared regulator po- Ta{1)) is the target nucleus thickness 
function. The exponential factor in Eq. (3) represents the probability that the 
parton suffered no semihard scatterings after the n-th one. In such a way, unitarity 
is explicitly implemented at the nuclear level, as discussed in Ref. [1,3]. The sum 
over n may be explicitly performed in Fourier space. The result reads: 

d'^PT 

where 



/ ^"^'''^''^^harArrlPo) , (4) 



SLrdirT;P0) = J d% [e-~<a..irT.P0)TAb) _ e-.L..(Pc)T.(6)] 



and 



<arArT;Po)= J d'k[l-e-^^'-~-]^ . (6) 

Note that <^^^^(rT) cx 4 as rr and crLrd(^T) ^ ^Lrd as rr ^ oo. This 
suggests the interpretation of <yfi„,.,i{r) as a dipole-nucleon "hard" cross section. 
This dipole is of mathematical origin, and comes from the square of the scattering 
amplitude written in the Fourier variable Tt, which represents the transverse size 
of the dipole. Then, we can interpret S^^^.^ as the dipole- nucleus "hard" cross 
section. Eq. (5) clearly incorporates Glauber-Gribov multiple scatterings of the 
colour dipole. Note that no nuclear effects on PDF's arc included, but shadowing 
is partly taken into account by the dipole multiple scatterings. 
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2. CGC and the colour dipole 

The Colour Glass Condensate is an effective theory for the nucleus gluon field at 
small-x, which describes the high-density regime where gluon saturation comes into 
play to modify the free-nucleon PDF's [7]. 

In this framework, minijet production in pA collisions may be related to the 
cross section for the scattering of a colour dipole on a nucleus [8]. The basic as- 
sumption is that the density of projectile partons is low enough for parton corre- 
lations inside the proton to be due only to DGLAP evoultion. In this case one is 
allowed to treat the proton in the coUinear factorization limit, exactly as in the 
pQCD-|-Glauber model of the previus section, see Eqs. (1) and (4) compared to 
Eqs. (36) and (32) of Ref. [8]. 

However, the dipole-miclcus interaction is computed entirely in the CGC model 
for the target nucleus. By using the so-called "Gaussian approximation" [9], in 
which parton correlation are assumed only to be Gaussian, one may write the dipole- 
nucleon cross section appearing in Eq. (4) as follows: 

al,,irr) = 4.a. N^J^^i^' e^'^-^^] , (7) 

where Nc is the number of colours, and t = ln(l/a;) is called "rapidity". Next,, fXr 
is interpreted as the unintegrated gluon distribution and has two limits [9] : 



47r2 dxG{x,k^) 
iV2-l dln(fc2) 
^ ,2, Q" 
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(8) 



These high- and Iow-Zct limits are defined relatively to the "saturation momentum" 
Q1{t) = Qge'*-^^"''^'^"'^"^ where tq is the minimum rapidity at which gluon sat- 
uration occurs and Qo is the corresponding saturation scale. Unfortunately, the 
constants 5, tq and Qq cannot be computed in this approximation. 



3. CGC is pQCD+Glauber 

The similarity of the CGC model at high-A:T with the pQCD-|-Glauber model be- 
comes evident by comparing Eqs. (7) and (8) with Eqs. (5), (6) and (2). To make 
this relationship more precise we proceed in four steps. 

Step 1. In the CGC, we separate "soft" and "hard" interactions: 

fj _ Att"^ dxG{x,k'^) k^ 



Air = /uf + /xf with \ d ln(yfc2 ) (k^ + Qlf (9) 

SO that fi^ vanishes quickly for kx < Qs- Accordingly, we have a soft and hard 
dipole- nucleon cross section: ^^q^ = 5^ + 5^ . 
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Step 2. In the dipole-nucleus cross section, we isolate contributions from pro- 
cesses with purely soft rescatterings, and processes with at leat one hard scattering. 
These define, respectively, the soft and hard dipole-nucleus cross section: 



The first term includes any number of soft scatterings but no hard ones, and the 
second term processes with at least one hard scattering. 

Step 3. Study observables for which it is possible to neglect the soft part, e.g., 
hadron spectra at large pr (Cronin effect), or integrated charged multiplicities at 
large ^/s, where contribution from hard processes should become dominant. 

Step 4. Following Ref. [10], we use the DGLAP equation to approximate the 
unintegrated gluon distribution in Eq. (9) with the integrated one: 



Eqs. (10) and (11) are equal to the pQCD-l-Glauber Eqs. (5) and (6) with po = Q, 



Under the assumptions stated in the previous section, it is possible to quantitatively 
compute some observables of interest. While a detailed applicaton of Eq. (4) to 
hadron p^-spectra and the Cronin effect is under investigation [11], we may use 
it to study RHIC Au-Au data on the centrality dependence of charged particle 
pseudorapidity densities [12]. The model was derived in the case oi pA collisions, 
but is generalizable to AB collisions under the assumption that in the regime of 
interest both nuclei are composed of a diluted enough system of partons. 

As a first step, we need the average minijet multiplicity at fixed impact pa- 
rameter and pseudorapidity. This is obtained by integrating Eq. (3) over pT, and 
introducing the thickness function Tb of projectile nucleus: 




In fcy TT 

This approximation is valid at large Q'^))Ql. 
Result. Collecting all the steps, we obtain 





and X — x^^^. 



4. What do RHIC data have to say? 




(po)Ta(0) +^^5 , 



(12) 



Saturation as Glauber-Gribov multiscatterings 



5 



S2 
T3 




400 



pQCD+Glauber 



PHOBOS 



Fig. 1 Charged particle mul- 
tiplicity per participant pair at 
RHIC, Eq. (13), as a function 
of the number of participants at 
different center of mass energies. 
Solid lines are the model results for 
the hard plus soft component. In 
the numerical computations both 
the gluons and the quarks have 
been included. Dashed lines are 
the result of the KLN saturation 
model [14]. Experimental data are 
from the PHOBOS collaboration 
[15]. 
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where the K-factor that simulates higher-order corrections to the pQCD cross sec- 
tion is exphcitly shown. 

The integrand in Eq. (12) may be interpreted as the average density of projec- 
tile partons (at a given x) times the probability of having at least one semi-hard 
scattering against the target. At low values of po the semihard cross section is large 
and the target becomes mor and more black to the projectile partons: the probabil- 
ity of scattering at least once becomes so high that nearly every projectile parton 
scatters and is extracted from the incoming nuclear wave-function. In this regime 
even if we use a lower cutoff po no more partons are there to be extracted. For this 
reason the minijet multiplicity tends to saturate [4]. We call saturation cutoff Psat 
the largest value of po at which this happens, and evaluate minijet multiplicities 
from Eq. (12) with po = Psat- We refer to [5, 12] for more details. 

To apply this minijet-level computation to integrated charged multiplicities, we 
write the charged particle multiplicity per unit rapidity as the sum of a soft and a 
semi-hard part [12]. For the soft part we use the wounded-nucleon model, which 
postulates a scaling of iV^'* with the number of participants. The semi-hard part 
is assumed to be completely computable from the saturation criterion for minijet 
production described above. To convert the minijet multiplicity to charged particle 
multiplicity, we further assume isentropic expansion of the initially produced minijet 
plasma, and parton-hadron duality. Our final formula reads 



1 



Npart{b)/2 dT] 



-{b) 



1 



5 d7V™J 



T-soft 



Npartib)/2 3 d?7 



ib]Psat) ■ 



(13) 



As the soft component n^o/t = nsoft{\fs) in Eq. (13) does not depend on the 
centrality of the collision, we may fit it to data for central collision only (we used data 
at Vs = 56-200 GeV [13] and extrapolated the fit down to 19 GeV, sec Rcf. [5]). The 
behaviour of the observable in non-central collisions is then completely determined 
by the model. In Fig. 1 we show the results of Eq. (13) compared to PHOBOS data. 
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Note that in our computation we used standard DGLAP evolved PDF's from the 
GRV group [6], with no nuclear modification of any sort. 

5. Conclusions 

The pQCD+Glauber formulation gives a very well defined framework in which de- 
tailed numerical calculations may be performed, and which is well-tested in the 
case of micleon-miclcon collisions. It is equivalent to the CGC model in a kinematic 
region where parton densities are low enough for gluon saturation effects to be 
negligible, but high enough for multiple scatterings and unitarity effects to become 
important. RHIC data on charged multiplicities at central rapidity are very satisfac- 
torily described by the pQCD+Glauber model with free-nucleon PDF's. Therefore, 
at least at central rapidity, RHIC data do not show any sign of gluon saturation. 

However, integrated multiplicities have very limited sensitivity to the details 
of the production processes, and alone do not allow any definitive conclusion. To 
assess the presence or absence of nuclear effects beyond multiscatterings and unitar- 
ity, it is necessary to study more differential observables in a cleaner environment 
than the hot and dense medium produced in Au-Au collisions. Novel effects may 
then be observed as a deviation of experimental data from the baseline given by 
pQCD+Glauber. The Cronin effect on hadron production in dA collisions at RHIC 
is an ideal candidate. 
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